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The CDF Collaboration has analyzed 955 pb~^ of CDF II data collected between March 2002 
and February 2006 to search for electroweak single top quark production at the Tevatron. 
We employ three different analysis techniques to search for a single top signal; multivariate 
likelihood functions; neural networks; the matrix element analysis technique. The sensitivities 
to a single top signal at the rate predicted by the Standard Model are 2.1 a, 2.6 a and 2.5 a, 
respectively. The first two analyses observe a deficit of single top-like events and set upper 
limits on the production cross section. The matrix element analysis observes a 2.3 o" single top 
excess and measures a combined t-channel and s-channel cross section of 2.7^13 P^- Using 
the same dataset, we have searched for non-Standard Model production of single top quarks 
through a heavy W' boson resonance. No evidence for a signal is observed. We exclude at the 
95 % C.L. W' boson production with masses of 760 GeV/c^ (790 GeV/c^) in case the right 
handed neutrino is smaller (larger) than the mass of the W' boson. 



1 Introduction 

In 1.96 TeV proton anti-proton collisions at the Tevatron, top quarks are predominantly pro- 
duced in pairs via the strong force. In addition, the Standard Model predicts single top quarks 
to be produced through an electroweak t- and s-channel exchange of a virtual W boson (Figure 
[T|). The production cross sections have been calculated at Next-to-Leading-Order (NLO). For a 
top quark mass of 175 GeV/c^ the results are 1.98±0.25 pb and 0.88±0.11 pb for the t-channel 
and s-channel process respectively^. The combined cross section is about 40% of the top anti- 
top pair production cross section. The precise measurement of the production cross section 
allows the direct extraction of the Cabibbo-Kobayashi-Maskawa matrix element \Vtb\ and offers 
a source of almost 100% polarized top quarks^. Moreover, the search for single top also probes 
exotic models beyond the Standard Model. New physics, like flavor-changing neutral currents 
or heavy W' bosons, could alter the observed production rate^. Finally, single top processes 



Figure 1: Leading order Feynman diagrams for s-channel (left) and t-channel (right) single top quark production. 



result in the same final state as the Standard Model Higgs boson process WH — > Wbb, which 
is one of the most promising low mass Higgs search channels at the Tevatron^!. Essentially, all 
analysis tools developed for the single top search can be used for this Higgs search. 

2 Event Selection 

Our single top event selection exploits the kinematic features of the signal final state, which 
contains a real W boson, one or two bottom quarks, and possibly additional jets. To reduce 
multi-jet backgrounds, the W originating from the top quark decay is required to have decayed 
leptonically. We demand therefore a high-energy electron or muon (£'j'(e) > 20 GeV, or Pt(/^) > 
20 GeV/c) and large missing transverse energy (MET) from the undetected neutrino MET>25 
GeV. Electrons are measured in the central and in the forward calorimeter, |?7| < 2.0. Exactly 
two jets with Et > 15 GeV and \r]\ < 2.8 are required to be present in the event. A large 
fraction of the backgrounds is removed by demanding at least one of these two jets to be tagged 
as a 6-quark jet by using displaced vertex information from the silicon vertex detector. The 
secondary vertex tagging algorithm identifies tracks associated with the jet originating from a 
vertex displaced from the primary vertex indicative of decay particles from relatively long lived 
B mesons. The backgrounds surviving these selections are tt, W + heavy-flavor jets, i.e. W + bb, 
W + cc, W + c and diboson events WW , WZ, and ZZ. Instrumental backgrounds originate 
from mis-tagged W + jets events {W events with light-flavor jets, i.e. with u, d, s-quark and 
gluon content, misidentified as heavy-flavor jets) and from non-T^ + jets events (multi-jet events 
where one jet is erroneously identified lepton). 

3 Background Estimate 

Estimating the background contribution after applying the event selection to the single top 
candidate sample is an elaborate process. NLO cross section calculations exist for diboson and 
tt production, thereby making the estimation of their contribution a relatively straightforward 
process. The main background contributions are from W + bb, W + cc and W + c + jets, as 
well as mis-tagged W + light quark jets. We determine the W + jets normalization from the 
data and estimate the fraction of the candidate events with heavy-flavor jets using ALPGEN 
Monte Carlo samples^, which were calibrated against multi-jet data^. The probability that a 
W -|- light-flavor jet is mis-tagged is parameterized using large statistics generic multi-jet data. 
The instrumental background contribution from non-VF events is estimated using side-band data 
with low missing transverse energy, devoid of any signal, and we subsequently extrapolate the 
contribution into the signal region with large missing transverse energy. The expected signal 
and background yield in the W + 2 jet sample is shown in Tabled] and graphically as a function 
of -|- jet multiplicity next to the table. Our analysis is performed in the W + 2 jet candidate 
sample with at least one tagged 6-jet. This data sample features the highest signal purity. Table 
[H however, demonstrates that the expected amount of single top events is much less compared 
to the large amount of expected backgrounds. In fact, the uncertainty on the backgrounds is 
larger than the expected signal, which renders a simple counting experiment impossible. This 
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Number of Events 


s-channel 


15.4±2.2 


t-channel 


22.4±3.6 


W + bb 


170.9±50.7 


W + CC 


63.5±19.9 


W + cj 


68.6±19.0 


Mistags 


136.1±19.7 


non-iy 


26.2±15.9 


Diboson 


13.7±1.9 


Z + jets 


11.9±4.4 


ti 


58.4±13.5 


Total prediction 


587.1±96.6 


Observed in data 
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Table 1: (Left) Expected signal and background yield in the W + 2 jet sample where at least on jet is tagged as 
a &-jet. (Right) Graphical sample composition as a function of the W + jets multiplicity. 



is exactly the reason why the search for single top quark production requires the best possible 
discrimination between signal and background processes and motivates the use of multivariate 
analysis tools. 



3.1 Neural Network Jet-Flavor Separation 

Mistags and W + charm events are a large class of background where no real 6-quark is present 
and amount to about 50% of the W + 2 jets data sample even after imposing the requirement 
that one jet is identified by the secondary vertex 6-tagger. We use a neural network tool which 
uses secondary vertex tracking information to distinguish 6-jet events from charm and light- 
flavor jet events. Figure [2] shows the distribution of this jet-flavor separating neural network for 
the 644 W + 2 jets candidate events. All three single top analyses use this neural network tool 
to improve the sensitivity of the analyses. 



4 Analysis Techniques 

No single kinematic distribution encodes all conceivable signal background separation. We 
use several sophisticated analysis techniques to combine information into a single discriminant 
distribution which is used to extract the single top content in data. 



4-1 Signal Significance and Discovery Potential 

To quantify the significance of a potential single top signal we use the CLs/CLb method de- 
veloped at LEpI3. In this approach, pseudo-experiments are generated from background only 
events (without single top) and from signal plus background events. We calculate the probabil- 
ity (p-value) of the background only pseudo-experiments to fluctuate to the observed result in 
data. A-priori we quote the expected sensitivity to a single top signal as the median p-value 
obtained from the signal -|- background pseudo-experiments. I.e. the quantity represents 50% 
of the generated signal plus background pseudo-experiments that had a p-value equal or greater 
than the expected p-value. All sources of systematic uncertainty are included in our statistical 
treatment and we consider correlation between normalization and discriminant shape changes 
due to sources of systematic uncertainty (e.g. the jet-energy-scale uncertainty). 
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Figure 2: (Left) Jet-flavor separating neural network applied to W + 2 jets data. (Right) Likelihood function 
discriminant for t-channel single top. The histogram is normalized to the expected event yield. 



4-2 Multivariate Likelihood Function Analysis 

The multivariate likelihood function analysis computes a joint probability that a given candi- 
date event originates from signal or background processes given a set of event characteristics 
xi, ■■jXn^^^. The likelihood ratio, as given in Equation [U is used to build a likelihood function 
discriminant for s-channel and t-channel single top. 

r(x, X ) = = '-^ dl 

llj=l Psig~rlU=l Pbkg ^^sig + ^^bkg 

The t-channel likelihood function discriminant is shown in Figure [2] which used seven input 
variables including the jet-flavor separating neural network described in the previous sub-section, 
the mass of the lepton, neutrino and the tagged jet Mii^b, a Matrix Element calculation from 
MadEvent^ and other event kinematic variables. The corresponding s-channel discriminant 
uses six input variables. The best expected p- value of 2.3% (2.1cr) is achieved by combining the 
t-channel and s-channel to a combined single top signal. The observed data show no indication 
of a single top signal and are compatible with a background-only hypothesis (p- value 58.5%). 
The upper limit on the combined single top cross section is 2.7 pb at the 95% C.L. The best fit 
for the combined s-channel and t-channel cross sections yields ag+t = O.S^q'^ pb. 



4-3 Neural Network Analysis 

The neural network analysis combines 26 kinematic or event shape related variables to a dis- 
criminant output between -1 for background-like events, and -|-1 for signal-like events. The five 
most important input variables are reported by the neural network package. They are: the 
output of the jet-flavor separator neural network, M/j^^, the dijet mass, the pseudo-rapidity of 
the untagged jet multiplied by the charge of the detected lepton and the multiplicity of soft 
jets in the event with 8 GeV< Et <15 GeV. To separate t- from s-channel single top quark 
production, two additional networks are trained and a simultaneous fit to both discriminants is 
performed. The combined search features an expected p- value of 0.56% (2. 60"). The observed 
p- value is 54.6%, providing no evidence for single top production. The corresponding upper 
limit on the cross section is 2.6 pb at the 95% C.L. The best fit yields as+t = O.Olto'o pb- 



4.4 Matrix Element Analysis 

Using the Matrix Element analysis technique, we compute event-by-event probability densities 
that a given candidate event resulted from a given underlying interaction (signal or background 
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Figure 3: (Left) Distribution of the Neural Network analysis discriminant compared to data in the signal region 
(Left). The histogram is normalized to the expected event yield. (Right) Event probability discriminant of the 
Matrix Element analysis. The histogram is normalized to the best fit value. The inset shows the signal region 

with high discriminant values. 



hypothesis). The measured four- vectors of the observed jets and the charged lepton serve as 
experimental input. The probabihty density is computed by integrating over the parton-level 
differential cross section du, which includes the leading order matrix element for the process 
(calculated using MadEvent^, the parton distribution functions f{xi), and the detector resolu- 
tions parameterized by transfer functions W{y, x). Lepton momenta and jet angles are assumed 
well measured while the jet energy measurements are corrected to parton level energies using 
jet-energy to parton-energy transfer functions. We integrate over the quark energies and over 
the z-momentum of the neutrino to create a final probability density. 

P{x) = da{y)dqidq2f{xi)f{x2)W{y,x) (2) 

We use these probability densities to construct a discriminant variable for each event (Equation 
E]). We also introduce extra non-kinematic information by using the output (b) of the neural 
network jet-flavor separator which assigns a probability (0 < 6 < 1) for each 6-tagged jet of 
originating from a b quark. 

]^PD — • Psingletop 

b ■ Psingletop + b ■ Pwbb + {I — b) ■ Pwcc + (1 - b) ■ Pwcj 

The expected p-value of the combined search is 0.6% (2.50"). The observed p-value is 1.0% 
{2.3a), providing a hint for a single top signal. The best fit yields cTs+t = '^■'^-15 pb- 

5 Search for Heavy W' Resonances 

Using the same dataset, we have searched for non-Standard Model production of single top 
quarks through a heavy W' boson resonance, pp — > W tb ^ Wjjthat appear in models with 
left-right symmetry, extra dimensions. Little Higgs, and topcolor^. We look for unexpected 
structure in the spectrum of the invariant mass of the reconstructed W boson and two leading 
jets {mwjj)- No evidence for resonant W' production is observed and we exclude at the 95 % 
C.L. a W with Standard Model coupling strength and masses of 760 GeV/c^ (790 GeV/c^) in 
case the right handed neutrino is smaller (larger) than the mass of the W boson. These new 
limits exceeds similar searches performed by CDF in Run I and DO in Run II of the Tevatron 
program^l^. 




Figure 4: Background model compared to the data for the mwjj distribution for W + 2 and 3 jet events (left). 
Expected and observed 95% C.L. limits on the W' boson mass assuming Standard Model coupling strength (right). 



6 Conclusions 

We have performed searches for electroweak single top quark production at the Tevatron using 
955 pb~^ of data collected with the CDF II detector. The sensitivities to a single top signal at 
the rate predicted by the Standard Model range between 2.1 .. 2.6 a and are the most sensitive 
to-date. The multivariate likelihood function and neural network analysis observe a deficit of 
single top- like events in the data. The matrix element analysis observes a 2.3 a single top 
excess consistent with the Standard Model expectation. Using pseudo-experiment techniques, 
we estimated the compatibility of the three analyses to about 1.2 % given the correlation of about 
60% and 70%. Extensive cross-checks have been performed to understand the different outcomes 
in data. At present, there is no evidence for the cause other than statistical fluctuations given 
that the analyses work in different ways and make different, analysis specific assumptions. The 
larger datasets of 2000 pb~^, already available to the CDF experiment, will clarify what the 
data are trying to tell us. 

Using the same dataset, we exclude at the 95 % C.L. non-Standard Model single top quark 
production through heavy W' resonances with masses of 760 GeV/c^ (790 GeV/c^) in case the 
right handed neutrino is smaller (larger) than the mass of the W boson. 
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